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A METHOD AND APPARATUS FOR PROVIDING REAL-TIME 
CALCULATION AND DISPLAY OF STRAIN IN ULTRASOUND IMAGING 

CROSS-REFERENCE TO RELATED APPLICATIONS 
Provisional Application No. 60 / 061 ,290. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 

DEVELOPMENT 

Not applicable. 

BACKGROUND OF THE INVENTION 

The present invention relates to diagnostic ultrasound systems which measure 
and image anatomical structures, and their movement. More particularly, the present 
invention relates to a signal processing method and apparatus for calculating and 
5 displaying of strain velocities, in localized parts of the image, to be used in ultrasonic 
imaging systems. 

Recently, within the field of ultrasound imaging, physicians have become 
interested in using tissue strain and strain velocity for clinical measurements. 

The term "strain" refers to a characteristic of material being examined. For 

10 example, the strain associated with muscle tissue corresponds to a ratio of the muscle 
tissues initial length and the change in muscle tissue length during a prescribed time 
interval. In ultrasound imaging, the rate of change of strain (e.g., strain rate, strain 
velocity, etc.) may be visually presented to a physician as a colorized 2-dimensional 
image, where variations in color correspond to different strain velocities. It has become 

15 apparent that the viability of a segment of the muscle is related to the amount of muscle 
strain and temporal behavior of the strain that is performed by, or is imposed on the 

1 
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muscle segment. Also, it has been determined that malignant tumors may be detected 
based on their resistance to compression. 

One application of real-time strain velocity imaging is in cardiology. The strain 
velocity gives a direct and quantitative measure for the ability of the myocardium to 
5 contract and relax. By imaging along the myocard from an apical view, the local strain 
velocity component along the long axis of the heart can be measured. Measuring the 
local strain velocity component gives information about the local shortening and 
lengthening of the heart wall. By imaging from the parasternal view, the strain velocity 
component perpendicular to the heart wall can be found. Finding the strain velocity 

10 component perpendicular to the heart wall gives information about the local thickening 
of the muscle. Wall thickening measured with M-mode or from the 2D image is a 
commonly used measure for muscle viability. With strain velocity imaging, a direct 
measure for this thickening is available. The strain velocity images can potentially add 
to the diagnosis of a number of cardiac disorders including, for example: 

15 Another application of strain velocity imaging is in heart transplants. Velocity 

variations inside the myocardium are important for the diagnosis of rejection after heart 
transplantation. The strain velocity images give a direct display of these velocity 
variations. 

Another application of strain velocity imaging is in noninvasive 
20 electrophysiology. The preferred embodiment describes techniques to image the local 
contraction/relaxation contributions with a high spatial and temporal resolution. Local 
contraction/relaxation information can be used to accurately determine the localization 
of, for example, where the mechanical movement in the heart chambers is activated 
based on a cross section just below the AV -plane. Furthermore, abberent conduction 
25 pathways (Wolf-Parkinson- White) from the atrium to the ventricle can be localized for 
later ablation. Even the depth inside myocard of these paths can be better localized 
with this invention in order to determine if the patient should be treated with catheter 
techniques or surgical techniques. 

Another application of strain velocity imaging is in measuring cardiac wall 
30 thickening. A well established methodology in cardiac diagnosis is to acquire a M- 
Mode image and to measure the wall thickening of myocardium during systole. The 
preferred embodiment provides techniques to take this wall thickening information and 
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measure it in real-time with a high precision in both the spatial and temporal domain. 
The high diagnostic relevance of the current wall thickening measurements indicates 
that the imaging modality described in this invention contains highly relevant 
information for cardiac diagnosis 
5 To understand strain velocity in more detail, it is assumed that an object of 

initial length Lq may be stretched or compressed or itself lengthens or contracts to a 
different length L. The one-dimensional strain, defined as 

represents a dimensionless description of the change. If the length L is considered to be 
10 a function of time, and is described as L(r) = ri(t)-r 2 (/), the temporal derivative of the 
strain, the strain velocity, can be found using the equation: 



£= ¥ (2) 

If the velocity, v of every point in the object is known, an equivalent definition of the 
15 strain velocity is: 

£ = — (3) 

dr 

These equations also provide a useful description of the deformation of the object. In 
Eq. 3, r is the spatial direction of the stretching or compression. The relation between 
Eq. 2 and Eq. 3 can be seen if the length L is defined as L(0 = r?(0- r i (0» ^ Lo=L(roX 

20 where r\ is the distance to one end of the object, and r 2 is the distance to the other. As 
illustrated in Eq. 3, the strain velocity is in fact the spatial gradient of the velocity. The 
strain velocity thus measures the rate of the deformation of the object If the strain 
velocity is zero, the shape of the object is not changing. If the strain velocity is 
positive, the length of the object is increasing, and if the strain velocity is negative, the 

25 length of the object is decreasing. Strain velocity is also known as rate-of- deformation, 
stretching, strain rate or velocity strain. 
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Strain imaging is currently an established research area in ultrasonic imaging. 
The degree of deformation in imaged structure may be estimated by correlation of 2D 
images obtained before and after a pressure increase. One disadvantage of estimating 
image deformation based on correlation of images is that the instantaneous value of the 
5 strain is not calculated nor displayed in real-time. The lack of a real-time capability is 
an important clinical disadvantage. For example, if strain imaging could be performed 
in real-time, strain imaging could be applied more effectively in cardiac ultrasound or 
could be used as an interactive inspection modality where anomalies in tissue 
compressibility can be visualized in real-time according to the pressure gradients that 

10 are applied to the imaged structures. 

A method of position tracking has been proposed to estimate the local 
myocardial strain velocity based on radio frequency (RF) M-Mode acquisitions. The 
position tracking method is described in H. Kanai, H. Hasegawa, N. Chubachi, Y. 
Koiwa, and M. Tanaka, "Noninvasive evaluation of local myocardial thickening and its 

15 color-coded imaging," IEEE Trans, on Ultrasonics, Ferroelectrics and Frequency 

Control, vol. 44, pp. 752-768, 1997. However, the method described in the Kanai et al. 
article has the disadvantages of poor temporal resolution and high computational cost, 
which render real-time imaging difficult and costly. Furthermore, the method 
described in the Kanai et al. article is a manual M-mode technique, not well suited to 

20 form the basis for real-time two-dimensional strain images. Also, the strain velocity is 
a derivative of a velocity estimate and is therefore very noise sensitive. The 
fundamental velocity aliasing problem that is inherent in tissue velocity imaging makes 
noise difficult to overcome because aliasing prevents the pulse repetition frequency 
from being set at a low enough rate to allow a large observation time. If the 

25 observation time could be increased, the noise robustness of the strain velocity images 
could be significantly improved. 

A need remains for an improved ultrasound system to overcome the above- 
identified difficulties. It is an object of the present invention to meet this need. 



WO 99/17660 PCT/IB98/01663 

SUMMARY OF THE INVENTION 

It is an object of the preferred embodiment of the present invention to provide 
an ultrasound system that is capable of calculating and displaying strain velocity. 

It is a further object of the preferred embodiment to provide an ultrasound 
5 system that is capable of calculating and displaying strain velocity in real time. 

It is a further object of the preferred embodiment to provide an ultrasound 
system that is capable of displaying strain velocity for a complete two-dimensional 
region. 

Il is another object of the preferred embodiment to provide an ultrasound system 
10 that displays strain velocity and yet overcomes the fundamental velocity aliasing 
problems that are inherent in tissue velocity imaging. 

It is another object of the preferred embodiment to provide an ultrasound system 
that provides real-time imaging of strain velocity with high noise robustness. 

It is another object of the preferred embodiment to provide an ultrasound system 
1 5 that provides real-time imaging of strain velocity with high temporal resolution. 

It is another object of the preferred embodiment to provide an ultrasound system 
that provides real-time imaging of strain velocity in a computationally efficient and 
economical manner. 

It is still another object of the preferred embodiment to provide an ultrasound 
20 system that displays strain velocity in a number of display formats. 

These and other objects of the present invention are provided by an ultrasound 
system and method for calculation and display of strain velocity in real time. Velocity 
of each point of the image is determined by measuring the pulse-to-pulse Doppler shift 
at sample points along a scanline and calculating tissue velocity based on the Doppler 
25 shift. Strain velocity is then estimated as the gradient of the tissue velocity along the 
corresponding scanline. The tissue velocity gradient is estimated by measuring the 
velocity difference over a predefined distance, or by linear regression methods to 
estimate the slope of the velocity variation along the scanline. 

A beam interleaving technique is used to lower the pulse repetition frequency 
30 (PRF) without reducing the total data acquisition time for the image. This allows 
accurate estimation of small velocity differences. However, a low PRF could cause 
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frequency aliasing when the Doppler shift exceeds PRF/2, and thus give large errors in 
the velocity estimate. 

One embodiment of the present invention provides a method to directly deduce 
the difference in Doppler shift between two points in the image, which is insensitive to 
5 the frequency aliasing effect. A significant improvement in the estimate of strain 

velocity is obtained without compromising a high frame-rate by deducing Doppler shift 
differences in a manner insensitive to frequency aliasing. 

The result of the estimated strain velocity images are presented in the same way 
as for tissue velocity imaging: as a color-coded video images, as color-coded M-mode 
10 or as time-variation curves for single sample volumes. The color-coded images may be 
a mixture of the B-mode tissue image with the strain velocity image. The B-mode 
tissue image provides an anatomical reference for the detected anomalies in the strain 
velocity image. 

Other objects, features, and advantages of the present invention will be apparent 
1 5 from the accompanying drawings and from the detailed description that follows below. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates an ultrasound system according to a preferred embodiment of 
the present invention. 

Fig. 2 illustrates exemplary signal packets xl and x2 from two range positions 
5 along a scanline according to a preferred embodiment of the present invention. 

Fig. 3 illustrates correlation values Rl and R2 in the complex plane, for two 
different pulse repetition frequencies (PRFs). 

Fig. 4 illustrates an example of an intravascular ultrasound image created by a 
preferred embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

A method and apparatus are described for generating strain velocity 
measurements in real time. In the following description, numerous specific details are 
set forth in order to provide a thorough understanding of the preferred embodiment of 
5 the present invention. It will be apparent, however, to one of ordinary skill in the art 
that the present invention may be practiced without these specific details. 

A block diagram for an ultrasound imaging system according to a preferred 
embodiment of the present invention is shown in Fig. 1 . An ultrasonic transducer 5 
emits a pulsed ultrasonic beam 6 into the body. A transmitter 3 drives the transducer 5 

10 to emit the pulsed ultrasonic beam 6. The ultrasonic pulses are backscattered from 
structures in the body, like muscular tissue, to produce echoes which return to and are 
detected by the transducer 5, see Fig. 1 . A receiver 2 detects the echoes. The echoes 
are passed from a receiver 2 to a complex demodulation stage 12 and a tissue 
processing stage 1 1 . The complex demodulation stage 1 2 demodulates the echo signals 

15 to form I, Q data pairs representative of echo signals. 

The demodulated I, Q data pairs are complex Doppler signals that are passed to 
the strain calculation stage 14, which carries out strain velocity calculations as 
explained below. The complex Doppler signal is associated with a selected range in the 
area of interest of the color image. A complex Doppler signal typically comprises a 

20 segment of data samples which is used to estimate the Doppler shift. The echo signals 
are also passed to the tissue processing stage 1 1, which performs processing such as B- 
mode processing to form a 2D or 3D image of the anatomical structure scanned. 

Next, the discussion turns to the techniques carried out by the strain calculation 
stage 14 to calculate strain velocity. The strain calculation stage 14 calculates the 

25 velocity SV(r) from an estimate of the tissue velocity v(r) in two points r and rf dr 
according to the following equation: 

SJ/(r) = (v(r + rfr)-v(r) ) / dr (4) 

The strain calculation stage 14 estimates the tissue velocity v(r) and v(rt-dr) 
based on the demodulated complex Doppler signals for echo signals associated with the 
30 data points r and r+dr. The strain calculation stage 14 may utilize any of several known 
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techniques to calculate the tissue velocities v(r) and v(r+dr). Alternatively, the strain 
calculation stage 14 may calculate the strain velocity SV(r), also referred to as the 
velocity gradient, based on a linear regression analysis of the estimated tissue 
velocities. 

5 Alternatively, the strain calculation stage 14 calculates the strain velocity 

directly, without first calculating the tissue velocities. The strain velocity may be 
calculated directly based on the phase-variation of the complex correlation function, 
which is used in the "autocorrelation method", described below. 

When calculating the strain velocity, the strain calculation stage 1 4 may utilize 

10 filtering techniques to obtain more robust estimates of the radial tissue velocity 
gradient. An example of such a filter is given in P. Saint-Marc, J. Chen, and G. 
Medioni. "Adaptive smoothing: A general tool for early vision " IEEE Trans, on 
Pattern Anal and Machine Intel!., 13(6), June 1991. The adaptive smoothing filters 
avoid erroneous velocity gradients between different structures like blood/tissue. The 

15 adaptive smoothing filters of the strain calculation stage 14 may also produce spatially 
piecewise constant strain velocity values. These properties of adaptive smoothing 
make the technique applicable for extensive spatial and temporal filtering in order to 
increase the robustness of the strain velocity value calculated by the strain calculation 
stage 14. 

20 Next, an embodiment is described whereby the system utilizes autocorrelation 

to estimate tissue velocities. Optionally, the correlation function may be averaged in 
space with neighboring ranges and/or vectors, to reduce estimation variance. The 
Doppler shift is then calculated from the angle of the averaged correlation function. 
The correlation function used to calculate Doppler shift may be calculated by 

25 any of several known functions. For instance, an autocorrelation method for ultrasound 
blood flow velocity measurements can be found in C. Kasai, K. Namekawa, A. 
Koyano. and R. Omoto, "Real-Time Two-Dimensional Blood Flow Imaging Using an 
Autocorrelation Technique" IEEE Trans. Sonics Ultras., vol. SU-32. pp. 458-464 
(1985). Another autocorrelation method for tissue velocity imaging is described in 

30 W.N. McDicken, G.R. Sutherland, and CM. Moran, "Color Doppler velocity imaging 
of the myocardium;' Ultrasound Med. Biol., vol. 18, nos 6/7, pp. 651-654, (1992). 

9 
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In more detail, the velocity v can be calculated by first obtaining the pulse-to- 
pulse correlation function from the complex Doppler signal packet x(r,t); t=l,..,N. 

P(r) = I.conj(x(r,t))x(r,t); 

R(r) = lconi(x(r,l))x(r,t + \y, (5) 

5 where conj() means complex conjugate. In equation (5) the formula for the Doppler 
signal power P, as well as the complex autocorrelation R are given. 

Optionally, before performing the correlation operation, the strain calculation 
stage 14 may subject the complex Doppler signal to clutter suppression, such as a 
regression wall filter. The purpose of the wall filter is to suppress signals from 

10 acoustical noise, which often have a lower Doppler shift than the signals from the 
measured object. It is preferable to avoid systematic errors in the velocity estimate, 
which may occur with filters of this form. These problems are discussed in H. Torp, 
"Clutter rejection filters in Color Flow Imaging: A theoretical Approach," IEEE Trans, 
on Ultrasound, Ferroelectrics, and Frequency control vo/. 44, 1997. The degree of 

15 clutter suppression may be adjustable, and in some situations, when the acoustical noise 
level is low, the clutter filter may be omitted. 

The velocity is obtained from the phase angle of R: 

v(r) = c 1 2fo phase(R(r)) 1 2nT (6) 

20 

where c=speed of sound, and T is the time between the consecutive pulses. Once the 
velocity v(r) is calculated, the foregoing correlation function is repeated for velocity 
v(r+dr), where dr is an incremental distance which is a parameter of the estimation 
algorithm. An increased value of dr may give a decreased spatial resolution, but 
25 improve the accuracy of the strain velocity estimate. Next, the strain velocity SV(r) is 
calculated according to equation (4) above. 
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It is preferable that, when using a correlation function, the velocity not exceed 
the Nyquist limit. This occurs when the velocity magnitude exceeds the Nyquist 
Velocity limit, corresponding to phase(R) = n 

5 Vnyqyuist = c/(4*fo*T) = c*PRF/(4fo) (7) 

Next, an alternative embodiment is described, in connection with Fig. 2 
whereby the strain calculation stage 14 directly measures the strain velocity. The strain 
velocity can be found directly from the phase difference between correlation values at 

10 two different points in the image separated by a distance dr, see Fig. 2. The correlation 
values Rl=R(r) and R2=R(r+dr) are calculated by the strain calculation stage 14 using 
equation (5). The receiver 2 records at least two continuous echo signals for the 
corresponding regions within the two range positions 21 and 22. The complex 
demodulator 12 converts the echo signals to signal packets 24 and 25 of data samples. 

15 The strain calculation stage 14 calculates complex correlation functions Rl and R2, and 
powers PI and P2 from the packets 24 and 25. The strain calculation stage 14 includes 
an arithmetic module 26 which calculates strain correlation S. Then the complex 
product between R2 and conjugate of Rl is formed. The resulting complex number S 
will then have a phase angle equal to the difference between the phase angle of R2 and 

20 Rl 

S(r) = conj(R(r))*R(rHdr) (8) 

The "Strain correlation function" S(r) may also be averaged both in the beam 
25 direction, and between neighboring beams, to reduce estimation variance. 

Fig. 2 shows a sector scan 20 containing two range positions 21 and 22 along a 
scanline 23. The "strain correlation" S is calculated as S=conj(Rl)*R2, the angle of S 
is proportional to the difference in Doppler shift between range positions 21 and 22. 
The strain calculation stage 14 also includes a strain velocity module 27 and a 
30 reliability index module 28. The strain calculation module 27 uses the following 
expression to calculate the strain velocity SV(r) from the strain correlation function: 
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SV(r) = c/(47tdrTfo) phase (S(r)) 



(9) 



10 



15 



20 



25 



Usually, the velocity difference v(r+dr) - v(r) is much lower than each 
individual velocity. Therefore, the Nyquist strain limit, which occurs when phase(S)=7i 



is much less likely attained. 

The advantage of using the direct method for strain velocity measurement is 
apparent in the following typical example for cardiac application: Maximum cardiac 
muscle velocity is typically 0.1 m/s, and max strain velocity is 2.0 (m/s)/m. With fo=4 
MHz, the PRF should be at least 1 kHz to avoid velocity aliasing. If dr=8 mm, 
minimum PRF to avoid strain velocity aliasing is 160 Hz. By lowering the PRF with a 
factor 6, a similar decrease in estimation error can be obtained. 

In Figs. 3A and 3B, the effect of decreasing the PRF is illustrated. Figs. 3 A and 
3B show the correlation values Rl and R2 in the complex plane, for two different 
PRF's. The circles indicate the random fluctuations of Rl and R2 around the true 
value. In Fig. 3B, the PRF is lowered by a factor 2. The angle of Rl and R2 will then 
increase by a factor 2, and aliasing error will occur (i.e., because phase(R) > rc). 
However, the angle difference will still give the correct Doppler shift difference (i.e., 
because phase(S) < 7i), and the random error in the SV estimate will decrease. 

In parts of the image, more chaotic movement may occur, which makes strain 
velocity difficult to measure. The reliability index module 28 corrects for chaotic 
movement by calculating an index indicating the reliability of the strain velocity 
calculation. 

The reliability index module may use the following reliability index ri(r): 



The reliability index is compared with a threshold, to remove noisy SV-values from the 
display. The reliability index may also be used to modulate (or modify) the color scale 
used to display the strain velocity, for instance, by decreasing the saturation of the color 
when the ri - value is low. 



SVnyquist= Vnyquist/dr 



(10) 



ri(r) = abs( S(r) ) / (P(r)*P(r+dr) ) 



(11) 
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The preferred embodiments of the present invention may be implemented in a 
system that employs a beam-interleaving. Velocity imaging requires two or more 
pulses in the same beam direction to be able to measure the velocity. The result is a 
packet of signal samples in each range point along the beam, and the packet size N 

5 equals the number of pulses. Thus, N>=2. This signal packet (after complex 

demodulation) is the Doppler signal. The time between the pulses is referred to as the 
pulse repetition time, and its inverse, the pulse repetition frequency (PRF). The total 
observation time for the moving objects is then N/PRF. Generally, the velocity 
estimation error decreases when the observation time increases. This can be achieved 

10 either by increasing the packet size, or decreasing the PRF. In both cases, the 

acquisition time, and the frame-rate is reduced. Beam interleaving is a technique where 
M pulses (interleave-size) are sent sequentially in different beam directions in N cycles. 
In this way, the PRF is lowered by a factor M, without increasing the total acquisition 
time for the image, and thus keeping the frame-rate constant, independent of PRF. 

1 5 The most extreme situation occurs when the interleave size M equals the 

number of beams in the image. In this case, the packet size can be set to N=l , and the 
Doppler signal can be obtained by picking one signal sample from each frame in the 
real-time sequence of frames. Then the PRF equals the frame rate, and a continuous 
stream of Doppler signals samples is available for analysis in each point in the image. 

20 The calculation of strain can then be performed on a "time-sliding" signal packet of 
arbitrary length, and with arbitrary degree of overlap between the signal packets. 

The preferred embodiments of the present invention may be implemented in a 
system that offers B-mode imaging. Separate pulses can be used to generate the B- 
mode tissue image to obtain a tissue image with higher resolution than the tissue 

25 velocity image or strain velocity image, i.e. short pulse and high number of scanlines. 

Alternatively, however, the B-mode image may be obtained from the same 
pulses as the strain velocity image. For low PRF requirement, the correlation can be 
performed frame-to-frame. The equivalent PRF will then be equal to the frame-rate. 
The frame-rate can be further increased by applying "multi-line acquisition" (MLA). 

30 This technique allows reception of more than one receiver beam for each transmitted 
pulse. Typically, a slightly broader ultrasonic beam is transmitted, and the receiver 
beam-former is set up to receive and separate the signals from 2 or more different beam 
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direction within the transmit beam opening angle. In this way, an increase in frame rate 
can be obtained by limiting the number of transmit beam directions. 

The preferred embodiments of the present invention may be implemented in a 
system that offers harmonic imaging. Harmonic imaging, (or octave imaging) where 

5 the second harmonic frequency of the transmitted signal frequency is used for display, 
has been shown to reduce noise artifacts in B-mode images. This technique can be 
combined with the preferred embodiments of the present invention in three ways: first 
by using the harmonic B-mode tissue image in combination with the strain image as 
described above, second, by using the harmonic part of the received signal as a basis 

10 for the strain velocity calculation or third, by combining these first and second ways. 

Intravascular imaging has so far been mainly limited to tissue imaging. The 
application of tissue velocity imaging in intravascular imaging can provide information 
about the movement of the vessel wall relative to a catheter. 

Fig. 4 shows an example of an intravascular ultrasound image. The catheter 40 

15 usually acquires a circular image 41 . Inside the image there will be an interface 42 

between the vessel lumen and the vessel wall. The sample points 43 and 44 correspond 
with points 21 and 22 in Fig. 2. Similarly, the sample points 45 and 46 correspond with 
points 21 and 22 in Fig. 2. The radial motion of a given sample point like 43 relative to 
the catheter 40 can be accurately measured by tissue velocity imaging. Furthermore, 

20 the location of the sample point 43 can be tracked radially as a function of time if the 
sample volume is repositioned according to the measured velocity multiplied with the 
time delay between consecutive image frames. If this is repeated for two points 
diametrically positioned on the vessel wall (43 and 45), one can make very sensitive 
estimates of the variations in vessel diameter. Similarly, variations in the vessel area 

25 can be monitored if the above analysis is repeated for all radial directions. The 
variations in vessel diameter/area can be used alone or together with auxiliary 
information like pressure to estimate important physiological parameters of the vessel. 
The strain velocity imaging described in this patent can be used to image the 
compressibility of the vessel wall. This can potentially be very important in the 

30 differentiation between various types of soft and hard plaque. 

Strain velocity imaging can also be performed on soft parts of the body by 
applying an increasing pressure on the probe against the skin. The tissue will then be 
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compressed, but hard regions will be compressed less than soft regions. This difference 
will show up in the strain velocity image. This type of imaging can for instance be 
performed when looking for tumors in the breast, the prostate, the thyroid or the liver, 
since tumors often are less compressible than normal tissue. 

In the foregoing specification the invention has been described with 
reference to specific exemplary embodiments thereof. It will, however, be evident that 
various modifications and changes may be made thereto without departing from the 
broader spirit and scope of the invention as set forth in the appended claims. The 
specification and drawings are, accordingly, to be regarding in an illustrative rather 
than restrictive sense. 
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CLAIMS: 

1 . A method for generating real-time strain velocity images comprising: 
acquiring echo signals (2) for a plurality of range positions (21, 22) along an 

ultrasonic beam (6) in an area of interest to cover a spatial region; 

estimating in real-time strain velocities (14) for said range positions (21,22) 
inside said spatial region; and 

displaying (7) estimated strain velocities for each range position at spatial 
coordinates on a display unit associated with said spatial region, to provide a live, real- 
time image of said strain velocities for said spatial region. 

2. The method according to claim 1 , wherein the step of estimating strain 
velocities (1 4) comprises: 

estimating tissue velocity for range positions (21, 22) along the ultrasonic beam 
(6) based on the echo signals; and 

calculating the strain velocity as a spatial derivative of the tissue velocity. 

3. The method according to claim 2 wherein the spatial derivative is found 
with a linear regression (14) of the tissue velocity for range positions (21, 22) along the 
ultrasonic beam (6). 

4. The method according to claim 1 5 wherein the step of estimating strain 
velocities (14) comprises: 

estimating tissue velocity for range positions (21, 22) along the ultrasonic beam 
(6) based on the echo signals; and 

calculating the strain velocity by determining a velocity difference between 
estimated tissue velocities associated with at least a first and second range positions 
(21, 22) and dividing the velocity difference by a distance between the first and second 
range positions (21 , 22). 

5. The method according to claim 1, wherein the step of estimating the 
strain velocity (14) comprises: 
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estimating a complex pulse-to-pulse correlation R(r) for a number of range 
positions (21, 22) along the ultrasonic beam (6) based on the echo signals; 

determining a strain correlation function, S(r), (26) over a radial distance dr 
according to an equation S(r)= conj(R(r))*R(r+-dr); and 

calculating the strain velocity (27) according to an equation SV(r)= c/(47idrTfo) 
phase (S(r)). 

6. The method according to claim 1, wherein the step of estimating the 
strain velocity (14) comprises: 

estimating a complex pulse-to-pulse correlation for a number of range positions 
(21 , 22) along the ultrasonic beam (6), based on the echo signals; 

calculating a strain correlation (26) function from at least two range positions 
separated by a given radia] distance; and 

calculating the strain velocity (27) based on the phase of the strain correlation 
function. 

7. The method according to claim 6 wherein the strain correlation function 
is given by multiplying the conjugate of the complex pulse-to-pulse correlation for a 
first range position by the complex pulse-to-pulse correlation for a second range 
position where said second range position is located the given radial distance from said 
first range position (26). 

8. The method according to claim 7 wherein the strain velocity is given by 
dividing a numerator defined as the product of the phase angle of the strain correlation 
function and the speed of sound by a denominator defined as the product of 4, n 9 the 
given radial distance, the ultrasound frequency and the time between consecutive pulses 
of said multiple of pulses (27). 
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9. The method according to claim 6 further comprising the step of 
temporally averaging the strain correlation function before calculating the strain 
velocity (14). 

1 0. The method according to claim 6 further comprising the step of 
spatially averaging the strain correlation function before calculating the strain velocity 
(14). 

1 1 . The method according to claim 1 wherein the step of displaying (7) the 
estimated strain velocity (14) comprises combining a B-mode tissue image (11) and the 
strain velocities at corresponding spatial coordinates in a color encoding. 

12. The method according to claim 1 wherein the step of displaying (7) the 
estimated strain velocity (14) comprises combining the M-mode tissue image (11) and 
the strain velocities at corresponding spatial coordinates in a color encoding. 

13. The method according to claims 1 wherein the step of displaying the 
estimated strain velocities includes displaying the estimated strain velocity (14) for a 
complete two-dimensional region (7). 

14. The method according to claim 1 wherein the estimated strain velocities 
are estimated for a three-dimensional region (14) and displayed (7) with a three- 
dimensional visualization technique. 

15. The method according to claim 1 wherein said spatial region is limited 
to a single point and the strain velocities are displayed (7) as a spectral or strain 
velocity curve versus time. 

16. The method according to claim 1 wherein a said acquiring echo signals 
is performed with a low pulse repetition frequency. 
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1 7. The method according to claim 1 6 wherein said low pulse repetition 
frequency is combined with a beam interleaving technique for increasing the frame-rate 
by acquiring information from different beam directions in the time interval between 
two consecutive pulses in a given direction. 

5 

1 8. The method according to claim 1 7 wherein multiple line acquisition 
techniques are applied to further increase the resulting frame-rate. 

19. The method according to claim 1 further comprising the steps of: 

10 calculating a strain reliability index (28) as a function of signal power for each 

sample point inside the said spatial region; and 

modifying the display (7) of the strain velocity according to the strain reliability 

index. 



1 5 20. The method according to claim 1 9 wherein the strain reliability index 

(28) given by dividing the absolute value of the product (26) of the conjugate of the 
complex pulse-to-pulse correlation for a first range position (21) and the complex 
pulse-to-pulse correlation for a second range position(22) by the signal power at one of 
said first or second range positions (21, 22). 

20 

2 1 . The method according to claim 1 wherein a harmonic of the echo signals 
is used to estimate the strain velocities (14). 



22. The method according to claim 1 wherein the strain velocity images are 
25 generated simultaneously with the application of an external means to produce a 

pressure gradient in the spatial region. 

23. The method according to claim 1 wherein the strain velocity images are 
generated while simultaneously applying changing pressure against the imaged spatial 

30 region. 



19 



WO 99/17660 PCTYIB98/01663 

24. The method according to claim 1 wherein the echo signals are acquired 
using intravascular ultrasound techniques (40) to assess regional changes in the vessel 
wall properties. 

5 25. In an ultrasound imaging system, a subsystem for measuring and 

displaying a rate of change of strain in tissue area of interest, comprising: 

a receiver (2) receiving echo signals corresponding to an area of interest; 
a demodulator (12) converting said echo signals to Doppler signals associated 
with range positions (21, 22) within an area of interest; 
10 a strain rate module (14) determining in real-time strain velocity for said range 

positions (21, 22) in the area of interest based on said Doppler signals; and 
a display displaying a strain velocity image for said range positions. 

26. An ultrasound imaging system according to claim 25 wherein the strain 
1 5 rate module ( 1 4) comprises : 

a tissue velocity estimation module estimating tissue velocity for range 
positions (21, 22) along an ultrasonic beam (6) based on the echo signals; and 

a strain velocity calculation module calculating the strain velocity as a spatial 
derivative of the tissue velocity. 

20 

27. An ultrasound imaging system according to claim 26 wherein the strain 
velocity calculation module comprises a linear regression unit calculating the spatial 
derivative by a linear regression of the tissue velocity for range positions (21, 22) along 
the ultrasonic beam (6). 

25 

28. An ultrasound imaging system according to claim 25 wherein the strain 
rate module (14) comprises: 

a pulse-to-pulse correlation module estimating a complex pulse-to-pulse 
correlation R(r) for a number of range positions (21, 22) along the ultrasonic beam (6) 
30 based on the echo signals; 

a strain correlation module (26) determining a strain correlation function, S(r), 
over a radial distance dr according to an equation S(r)= conjR(r)*R(r+dr); and 
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a strain velocity calculation module (27) calculating the strain velocity 
according to an equation SV(r)= c/(47idrTfo) phase (S(r)). 

29. An ultrasound imaging system according to claim 25 wherein the strain 
5 rate module (14) comprises: 

a pulse-to-pulse correlation module estimating a complex pulse-to-pulse 
correlation for a number of range positions (21,22) along the ultrasonic beam (6), based 
on the echo signals; 

a strain correlation module (26) calculating a strain correlation function from at 
10 least two range positions (21, 22) separated by a given radial distance; and 

a strain velocity calculation module (27) calculating the strain velocity based on 
the phase of the strain correlation function. 



30. An ultrasound imaging system according to claim 29 wherein the strain 
15 rate module (14) further comprises a temporal averaging module temporally averaging 
the strain correlation function before calculating the strain velocity. 



31. An ultrasound imaging system according to claim 29 wherein the strain 
rate module (14) further comprises a spatial averaging module spatially averaging the 

20 strain correlation function before calculating the strain velocity. 

32. An ultrasound imaging system according to claim 25 further comprising: 
a strain reliability module (28) calculating a strain reliability index as a function 

of signal power for each sample point inside the said spatial region; and 
25 a display modification module (7) modifying the display of the strain velocity 

according to the strain reliability index. 
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